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Abstract 

We study the renormalization of normal mixing matrices, which includes hermitian and 
unitary matrices as particular cases. We give a minimal, multiplicative parametrization of 
counterterms, and compute the renormalized Lagrangian to one-loop order in several simple 
models with N species of fermions, both in on-shcll and MS schemes. In on-shell scheme the 



mass-degenerate case is considered separately. 

1 Introduction 



In theories with many particle species which mix non-trivially due to interactions, the mixing matrix 
generally requires renormalization like any other parameters in the Lagrangian. Such renormaliza- 
tion has been considered in pQ within the context of the Standard Model (SM), and in extensions 
Qh| of the SM with Majorana neutrinos in [2J. 

In this paper we consider mixing-matrix renormalization in a more generic setting. Specif- 
ically, we study the renormalization of normal mixing matrices (i.e., matrices commuting with 
their adjoint), which includes hermitian and unitary matrices as particular cases. We give a min- 
imal, multiplicative parametrization of counterterms, and compute the renormalized Lagrangian 
^ ■ to one-loop order in several simple models with N species of fermions, both in on-shell (OS) and 

MS schemes. In on-shell scheme the mass-degenerate case is considered separately. We work in 
dimensional regularization [3] throughout this paper. 

Mixing-matrix renormalization in the SM, and in its extensions and low-energy effective the- 
ories, is usually closely related to other issues such as renormalization of theories with unstable 
particles, gauge invariance of the renormalization procedure, and CP violation. The latter is out 
of the scope of this paper. We assume all particles to be stable. If that were not the case, one-loop 
self-energy parts should be replaced by their dispersive parts in OS renormalization conditions [3]. 
Gauge invariance of the renormalized CKM matrix in the SM has been considered in [E1E]- In our 
case, Abelian gauge invariance plays a role in the discussion of unitary mixing matrices in 

In the next section we consider a model with N fermion flavors coupled to a scalar particle 
through a Yukawa interaction, specified by a hermitian coupling matrix. Renormalization of this 
model, and in particular the structure of its counterterms, is considered in detail and computed 
at the one-loop level. The results are generalized to normal matrices in The particular case 
of unitary interaction matrices is treated in £j3J where the relation between our parametrization 
of counterterms and the one commonly in use in the literature is discussed. In JDwe give some 
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final remarks. We gather material relevant to all sections in four appendices. In appendix EI 
in particular, we give a parametrization for mappings of normal matrices that we find useful in 
discussing renormalization of mass and mixing matrices. 

2 Hermitian mixing matrix 

The simplest model of fermion mixing consists of N Dirac fields ij} a , a = 1, . . . , N, which we gather 
together in a column field ij), interacting with a scalar field <p through a Yukawa coupling. The 
Lagrangian is given by, 

£ = -l0 o (D + m^ - M )V + ^ H ^J - |0 O 3 - ^ (1) 

the subindex indicating bare fields and parameters. The Yukawa interaction term is specified 
by a matrix of couplings H , whose elements are the expansion parameters in perturbation theory. 
M is the fermion mass matrix, which is assumed to be regular (i.e., no two of its eigenvalues are 
equal). The degenerate case will be considered below in EJ2.3I If [M, H] ^ 0, the interaction mixes 
flavors. Clearly, M and H must be hermitian. The interaction terms in 3 and 4 are needed for 
renormalizability. If fermions are massless and £ = 0, C is invariant under the discrete symmetry 
ip — > 75^, <p ~~ ~4 > which forbids a (f) 3 term. For massive fermions such a term cannot be avoided. 
Because of its simplicity and of its lack of symmetries preventing renormalization of its parameters, 
C is in some sense archetypal, so we consider its renormalization in detail below. 

We never use the summation convention for flavor indices a,b,c, . . . . Space-time indices are 
denoted by Greek letters, with summation over repeated indices always understood. 

We write C in terms of renormalized fields and couplings by introducing renormalization 

1/2 

constants. The field <fr renormalizes multiplicatively <j) = <j) = (1 + l/2 5Z ( f ) )4>. In MS 
Z^ = Z^H, A, d)[7j, whereas in OS Z^ = Z^(H, A, d; rrv^/fj, 2 , M/fJ-, £///)• We do not expect mul- 
tiplicative renormalization of A, which already at one loop receives 0(H 4 ) divergent contributions 
from a box diagram, A = fi e (X + <5A). The dependence of 5X on the parameters of the theory is 
analogous to that of Z§. In these expressions d = 4 — e is the dimension of space-time and /i the 
mass scale of dimensional regularization. 

The mass will mix under renormalization with the other dimensionful parameters in C. By 
dimensional analysis, 

m4 = Zg\ mj + Szg^MatMci + dZ^f. 

abed 

In OS renormalized masses are physical, so M a b = rria^ab and, 

ab 

(i) 

the renormalization constants SZm^ depending on H, A, d and ratios of dimensionful parameters. 

( 7 ) TT71 

In MS there is no dependence of SZ^ on masses and £ [7j, but M is not diagonal. We can choose 
our flavor basis, however, so that at tree level M is diagonal. Off-diagonal elements in M are 
therefore of second order in H, A. At one loop we can then write, 

5ml = 5Z &\ m l + E SZ ab m ° m t> + SZ £l? ( 2 ) 

ab 
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in both schemes, masses being physical in OS and renormalized ones in MS . Similarly, for the 
cubic coupling we have, 

So = S /2 ZPt + ^ 2 Y, SZ Z% M ^ ( 3 ) 
ab 

where we can set = for a ^ b at one loop. At that order, £ mixes with M through a triangle 
diagram that does not depend on £ or m|. 

The renormalized fermion field can be related to the bare one as ip = Arj), with A a complex 
N x N matrix. It is more convenient to introduce the polar decomposition of A explicitly, 

V = C/Z 1 /V, U = e~ iSU , Z 1 ' 2 = l + Uz. (4) 

Both SU, SZ are hermitian. It is not difficult to show (see appendix [BJ), however, that we can 
always parametrize a unitary matrix U in a neighborhood of the identity as U = e~ lSU ' e~ lSU , with 
SU' a linear combination of diagonal generators and SU of the remaining ones. The effect of SU' is 
a flavor-dependent phase reparametrization of the fermion fields. Although C is not invariant under 
such transformations, physical quantities (such as S"-matrix elements) remain unaffected by them. 
We choose to set SU' = 0. In fact, as we shall see below, in OS renormalization conditions do not 
determine the diagonal elements of SU. (Thus, in Q we have SZ G u(N), whereas SU E u(N)/a, 
with u(N) the Lie algebra of N x N hermitian matrices, and a its Cartan subalgebra of diagonal 
generators.) 

The renormalization matrices U and Z 1 / 2 have different roles in the theory. If, for instance, 
we take tp a to be scalar fields and set A = in £, Z 1 / 2 turns out to be finite, as we would expect 
of wave-function renormalization in a superrenormalizable theory 8 . On the other hand, in OS, 
mass renormalization contributions make U divergent. If, furthermore, the scalar fields are real, 
Z must be real symmetric and U orthogonal. In this case iSU is antisymmetric, with diagonal 
elements vanishing by definition, which is consistent with the fact that no phase redefinitions are 
possible for real fields. 

The fermion mass matrix mixes under renormalization with £. Following the parametrization 
given in appendix El we write, 

M = Ul(M + SM) U m + 6Z®t, \SM, M) =0, U m = e~ iSU ™, (5) 

with SM, SU m , SZ^} hermitian. Here we have conventionally chosen an additive parametrization 
for the first term in M . Since we treat £ as a perturbative parameter, we collect all dependence 
on it in the second term. At one loop SZ$ = 0, so we can write M = M + SM with SM = 
SM + i [SU m , M]. SU m is only determined up to addition of a matrix commuting with M. 

Finally, the matrix of Yukawa couplings H is multiplicatively renormalized. The most general 
multiplicative transformation that can be applied to H preserving hermitianity is of the form (see 
appendix 0) , 

H = if/ 2 W*Z H HW, [Z H ,H}=0, W = e~ iSW , (6) 

with Zh-i 8W hermitian. In MS Zu and W can depend only on H, A and d, so that W commutes 
with both Zh and H, and drops from eq. (JBJ). This would not be the case, though, if there were 
another dimensionless matrix in C on which W could depend such as, e.g., another mixing matrix 
for an additional family of fermions. 

Substituting these definitions into Q we see that, out of the three unitary matrices we have 
introduced, only two combinations enter C, namely, U m U and WU. We can therefore always 
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choose U = 1. In OS, however, we parametrize the theory so that M and M are simultaneously 
diagonal by setting U m = 1. 

The Lagrangian can then be written in terms of renormalized fields and counterterms as, 



£ = ~<j>(D + m\)4> ~ i^Z^Ucj) - ^Am 2 </> 2 
+ V>0' - M)ip + tpSZi @ij> - ipAMip 
+ /i e/2 ^ (H + AH ) 
- i/W 2 (£ + A^ 3 - ^(A + AA)<^ 4 , 



(7) 



with the one-loop counterterms defined as, 

Am 2 = SZ^m^ + Sm 2 ^ ; AA = 5X + 2A5Z0 (8a) 
A£ = (5£ + :W; ^ = ^ (1) C + ^^gm a (8b) 



AM = - {M, SZ} + SM + i [517 + SU m , M] (8c) 

AiJ = SZ H H + i[SU + SW, H] + ^ {SZ, if} + hsZ+H. (8d) 

In order to compute the renormalization parts of the theory we have first to fix a renormalization 
scheme. 



2.1 Overspecified counterterm parametrization 



In © and (jHJ) we have given the relation between bare and renormalized hermitian matrices 
using the parametrization in appendix which is applicable to generic normal matrices. In the 
particular case of hermitian matrices a more obvious multiplicative parametrization is a congruence 
transformation H = A^ HA, with A non-singular and dependent on H. Such parametrization 
must satisfy constraint relations, since it is overspecified in the following sense. 

Let / : u(N) — ► u(N) be a map of hermitian matrices, given by /(JET) = A^HA, with 
A = A(H) non-singular. Then, we can find a non-singular B = B(H) such that, f(H) = B^HB 
and, 



bbKh 



0, 



0. 



(9) 



To see this, we notice that according to appendix there must be a hermitian matrix Zh com- 
muting with H and a unitary matrix U such that f(H) = ZfjHU. Defining B = Z l Jj 2 U , we 
obtain the result. The parametrization given in Appendix E] is minimal in the sense that eqs. Q 
are satisfied identically. 



2.2 On-shell scheme 

In OS we set U m = 1. If M is regular and diagonal, and [M,SM] = 0, then 5M must be 
diagonal, as well as M . 

The renormalization conditions we impose on the scalar field self-energy IL^, (p 2 ) are Il0(m 2 ) = 
= IL'.(mV), where the prime stands for d/dp 2 . Defining, 

ft<t>(p 2 ) = -g-f H ab H ba (m 2 a (m 2 ) + m 2 a (m 2 ) + ((m a + m b ) 2 - p 2 ) b (p 2 , m 2 , m 2 )) , (10) 

a, ft 
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which is the (D(e ) part of the unrenormalized (f) self-energy, we obtain, 

U^p 2 ) = ^(p 2 ) - n^m 2 ) - {p 2 - m 2 )n'4m 2 ) (11a) 

5Z + = "i Tr( ^ 2) + (llb) 
Sml = - JL- {2Tr(H 2 M 2 ) + Tr {{HM) 2 )) - f^(m 2 ) + m^m 2 ,). (11c) 

In (jl()j) . ao and 60 refer to finite parts of loop integrals, defined in appendix El (All dependence 
on the renormalization scale \i is through these integrals and, as is easy to check, li^p 2 ) does not 
depend on /x as it should in OS.) All masses in this section are physical (pole) masses. 
For the fermion two-point function, which is a matrix in flavor space, we write, 

r=yi-M -n(p), n(p) =^ v {p 2 ) + ^ s (p 2 ), (12a) 

where the form factors are given by, 

Ty(p 2 ) = --l-H 2 -8Z + n v {p 2 ), ^p 2 ) = --± r HMH + AM + n s (p 2 ), (12b) 
with, 

tta b (p 2 ) = -^-2 Yl H acH cb b„(p 2 ,ml, m 2 c ), £l S ab {p 2 ) = ^2 Yl H acH cb m c b (p 2 , mj, m 2 c ). (12c) 

c c 

Renormalization conditions are expressed in terms of 'S V ' S (p 2 ) as (I], 

m a Yy aa (m 2 a ) + Sf a (ml) = 0, £^(m 2 ) + 2m 2 £^(m 2 ) + 2m a Z*(m 2 a ) = (13a) 

m b ^ b (m 2 ) + £f 6 (m 2 ) = 0, ™a^(™a) + £f b (™a) = 0, (13b) 

where in the second line a ^ b. With these conditions, we obtain for the diagonal two-point 
functions, 

^L(P 2 ) = ^L(P 2 ) ~ ^LK) " Zrnln^ml) - 2m a n s a ' a (m 2 a ) (14a) 
£f> 2 ) = V s aa (p 2 ) - n s aa (ml) + 2mX(m 2 a ) + 2m 2 ^^(m 2 J, (14b) 

and for the off-diagonal ones, (a ^ b) 

ZUP 2 ) = ^(P 2 ) " m ^ K) " m ^ K) - ^IziSW (14c) 

m a -m b m a -m b 

^) = ni(p 2 ) + (n&<mS) - + "' t " ft(m ° ) - m °" ftK) . d4d) 

"ia — w-fe m a — m b 

In both cases, diagonal and off-diagonal, renormalization constants can be compactly expressed as, 

SZ = --^ r H 2 + ft v (p 2 )-^ v (p 2 ) ; AM = —^r-HMH — Q, s {p 2 ) + H s (p 2 ) . (15) 

It is immediate to check, by substituting the expressions for H v,s into these equations, that SZ 
and AM are hermitian. With this, the fermion propagator is renormalized to one loop. We can, 



5 



however, compute the mass counterterm SM and SU from the expressions Q15|) for AM and SZ. 
From (|5c|) we have, 

SM + i[SU,M} = AM -^{SZ,M}. (16) 

This equation can always be solved for SM and SU, as discussed in appendix[0 In this case, in 
which M is diagonal and regular, however, (|16|) is trivial to solve. SM is given by the diagonal 
elements of the rhs of (|16|) and SU by the off-diagonal ones, 

5M ab = (AMaa ~ m a SZ aa )5 ab (17a) 

SUat = *- ( AM ab - ^±J^L5Z ab \ ■ a^b. (17b) 

m a -m b \ 2 J 

The diagonal elements 5U aa are not determined by renormalization conditions. As mentioned 
above, they only change the phases of fermion fields, and can be chosen to vanish. 

We consider next the renormalization of the Yukawa coupling, given by the 1-PI three-point 
Green's function T(pi,p2). (Here p\ is the momentum of the incoming scalar and j>2,3 are the 
momenta of the outgoing fermions.) At tree level V = H is a Lorentz scalar. Expanding the 
one-loop r in the usual 7-matrix basis, it is clear that only the scalar form-factor can receive 
divergent contributions, since the counterterm for T in C is a scalar, the other form factors being 
finite. This is easily seen also from the explicit form of the corresponding loop integrals, since by 
power counting only those terms containing two powers of loop momentum in the numerator are 
divergent, and they contribute to the scalar form factor only. 

We will then focus on F(p\,p2) = l/4Tr(r(pi,2?2)). For concreteness, we assume that the 
physical value of the coupling H is fixed by this form factor. F is a scalar function of momenta, 
depending on pi, P2 only through p\ , p\, (pi — pi) 2 ■ These combinations are fixed if the external 
momenta are required to be on their mass shell. Calling F^ the value of F ab (p\,p2) at p\ = m\, 
p\ = m 2 , (pi — P2) 2 = in?, we impose the renormalization condition 

i> s ) = H. (18) 

(We notice, furthermore, that the hermitianity of the lhs is guaranteed by CPT invariance, which 
also requires H to be hermitian.) At one loop F(p\,p2) is given by 

F ab { Pl ,p 2 ) = fi e/2 (H ab + AH ab ) - fi e / 2 ^- (H 3 ) ab + fi e / 2 Y,H a cH cd H db h 1 (p 1 ,p 2 ;ml,m 2 c ,m 2 d ) 

c,d 

~ ^ e/2 ^^H ac H cb h 2 (pi,P2-,ml,rn 2 c ), (19) 

c 

where the loop integrals h\p are defined in appendix [DJ The renormalization condition 1)18(1 then 
leads to, 

AHab = 8^~e ( H3 ^b ~ 2 H acH cd H db h { ° s) (m 2 , m 2 c ,m 2 ) + £ £ H ac H cb h 2 ° s) (m 2 , m 2 c ), (20) 

cd c 

where lif^ refer to the on-shell values of those integrals. This counterterm is enough to renormalize 
the scalar form-factor at one loop. We can, in principle, determine the counterterms to the coupling 
matrix H as defined in (|8d|) through the equation, 

SZ H H + i [SW, H] = AH - i [SU, H] - ~ {SZ, if} - ^SZ^H, (21) 
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where all the quantities on the rhs are already known. Equation (|21|) always has a solution, since 
by definition [SZ h, H] = 0, as shown in appendix|0 Unlike (|16j) . in this case it is more difficult to 
find the solution algebraically in closed form, only some of the contributions to each counterterm 
being obvious, 

6Z H = JL (if 2 ) - IsZ* + • • • ; [SW,H] = - [SU,H] + ■■■ . (22) 

Eq. (j2"T]) can be solved by projecting it over an appropriate basis for the algebra u(N) (see appendix 
which can be done numerically. It should be clear, however, that counterterms are completely 
fixed by renormalization conditions. Once those are established, no other choices are involved. We 
notice that diagonal elements of SU contribute to SW, as expected, since a change of phase of 
fermion fields leads to a corresponding change in H. 

2.3 The mass-degenerate case in on-shell scheme 

As discussed in 32~2l in OS we set U m = 1, so that M = M + SM, with [M,SM] = 0. M 
and M can then be simultaneously diagonalized. If M is regular, choosing a basis in which it 
is diagonal immediately implies that 8M is diagonal, and then so is M Q . If M is degenerate, 
though, 8M need not be diagonal even if M is. We notice also that H can be diagonalized within 
each eigenspace of M by means of a unitary transformation of O(H ). We assume that such a 
transformation has been performed and then, in general, SM will not be diagonal within mass 
eigenspaces. 

For notational simplicity we assume that there is one degenerate mass eigenvalue, say mi, with 
multiplicity l<r<N(r = N being the trivial case M oc 1), all other masses, m r+ i, . . . , itin, 
being non-degenerate. The generalization to the case of several degenerate eigenvalues presents no 
difficulties. 

The renormalized fermion two-point Green's functions, as given by 1|14|) and ()15|) are unchanged 
in the mass-degenerate case if a = b or if a ^ b but max(a, b) > r. If a ^ b and both a,b < r, those 
results need modification. 

Renormalization conditions 1)1 3b|) reduce to the single equation (a ^ b, a,b < r) 

mi££(m?) + £fb(m?) = 0, (23) 

leading to, 

AM ab - mi 5Z ab = J^- {H 2 ) ab + ^-(HMH) ab - m^mj) - ftf fc (m 2 ). (24) 

From definition (|8c|l we see that in this case, AM a b — mi5Z a b = 8M ab . Thus, the renormalization 
condition (|2"3*)l fixes the off-diagonal counterterm 5M a b to the value given in (|24[). and that is all 
that is required to renormalize the two-point function within the eigenspace of mi. SM a i,, as given 
by (|17a|) when a = b or a ^ b, max(a, b) > r, and by (j24j) . when a ^ b and a,b < r, obviously 
commutes with M, since it is non-diagonal only within eigenspaces of M. Furthermore, we set 
5U a b = 0, for a,b < r, instead of (|17b|) . Notice that when r = N, H is diagonal and therefore so 
are ft V ' S and SM. 

V s 

If, however, we want to make each form factor T, ab ' finite separately, we may impose additional 
renormalization conditions. These are quite arbitrary, as long as they are consistent with (|23j) (or 
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(|24jl). We can, for instance, take the limit of degenerate masses in (|14j) and to get, for a ^ b, 
a,b < r, 

SZ ab = -^(J^U + VL v ah {m\) + 2m?n£(mi) + 2™i^(™?) (25a) 
AM a6 = ^-(HMHU - n s ab (ml) + 2m?^'(m 2 ) + 2m?fiS(m?). (25b) 
Requiring instead Y^ b s (m 2 ) = 0, leads to, 



SZ ab = --l^(H 2 ) ab + n v ab {m\) ; AM a6 = -^{HMH) ab - tt s ab {m 2 ). (26) 



Both ((251) and (J2EJ) are consistent with 
2.4 MS scheme 

In MS we set U = 1. We choose a flavor basis for fermion fields so that at tree level the renormalized 
mass matrix M is diagonal, M = _M" p h ys - Off-diagonal elements of M are then 0(H 2 ). We write 
M = M 1 + iW, with M' = diag(mi, . . . , mjv) containing the renormalized masses and M the 
off-diagonal elements, and treat M as an interaction term. We write the fermion Lagrangian as, 

C,p = ^(i$- M')il> - ^MV + ^(SZi $ - AM)^ 

AM = 8M + i[SU m , M] + -{M, 8Z}, 

instead of the second line of 0. The tree-level fermion propagator is thus flavor diagonal. 

We notice, parenthetically, that M can be written as M = exp(iE)M' exp(— iE), with M' 
the diagonal matrix of eigenvalues and exp(— iE) the unitary matrix of eigenvectors of M. Our 
choice of tree-level flavor basis implies E = 0(H 2 ). At one-loop level, then, M' is given by the 
diagonal entries of M and M = i[E, M'] contains the off-diagonal ones. 

Counterterms for the scalar two-point function, as defined by ® and (jHJ), are given by the 
e-pole terms of (jllb|) and Ullcj) . The 4> self-energy is then H^lp 2 ) = f^p 2 ), with defined 
in l)10|) . Requiring that the two-point function have a zero at p 2 = m^p^ leads to the relation 

m l = m l P h ~ n 4>( m lph) at one lo °P- 

The fermion two-point function can also be read off the corresponding OS results. Starting 

from CEJ), we set SZ and AM to be the e poles of ((EI), to obtain S 7 ' 5 ^ 2 ) = ft V ' S (p 2 ). The 

counterterms SM and SU m can then be obtained from SZ and SM in the same way as in OS 

(see (|16jl ). In this case, unlike in OS, they are not needed to renormalize the Yukawa coupling. We 

define 

SM = AM - i{M, SZ} = (hMH + i {M, H 2 }^j , (27) 

which we will use in the one-loop renormalization group (RG) equation for M below. 

We require r aa (p)ii a (p)L2_ m 2 — in order to obtain renormalized masses in terms of physical 

" aph 

ones, leading to the relation 



m n 



m aph - {m aph n^ a (m 2 aph ) + ^f a (m 2 ph )) . (28) 



Expressing the off-diagonal elements M ab , a ^ b, in terms of physical masses and coupling constants 
involves a choice of parametrization of the theory. In principle, any value of M consistent with the 
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renormalization group equations (given below) is admissible. Different choices of M at one-loop 
level will result in different parametrizations of the renormalized masses m a (the eigenvalues of 
M), in terms of the physical ones m ap h at two loops. We could give M a definite value at a mass 
scale /Z , and evolve it with the RG equations to the desired scale Jl. For concreteness, we quote 
the expression for M obtained by matching the theory in MS to the OS results, 

M ab = (M ph ) afe - - {m aph n^ b (m 2 aph ) + n s ab (m 2 aph ) + m bph n^ b (m 2 bph ) + n s ab (m 2 bph )) , (29) 



where M" p h = diag(mi p h, • • • , mjvph)- When a = b this equation reduces to (J2i 

It is convenient to discuss at this point the effect of a finite, unitary renormalization of fermion 
fields, with SU = 0(H 2 ) and ^-independent. (We exclude from consideration unitary transforma- 
tions with SU = 0(H ), which correspond to transformations of the classical fields.) The effect 
of such a transformation on the basic Lagrangian is to change M U ] MU and // U ] HU. 
It is clear that M remains diagonal at O(H ), so our choice of a tree- level flavor basis is not 
altered. Also, M' remains unchanged through 0(H 2 ), so (|28[) still holds. Changes in M' start 
at 0(H 4 ). On the other hand, M changes by a term i[6U,M'] at 0(H 2 ). Counterterms still 
are of the MS form after the transformation, which is multiplicative and independent of e. (If, 
however, SU depends on masses, counterterms acquire mass dependence.) We see, then, that this 
U(N) freedom to perform finite unitary renormalizations of 1/7 is a source of arbitrariness in M. 
By choosing a definite value for M, as in ()29j) . we are reducing the ambiguity in the choice of 
one-loop flavor basis for if) to the subgroup of U(N) that leaves M invariant. If M is regular, this 
is the {U{1)) N Abelian subgroup of flavor-dependent phase transformations. 

The relation between bare and renormalized Yukawa couplings is given by ©• As discussed in 
relation to that equation, the unitary renormalization matrix W is trivial in MS in this model. 
From the value of A.H in OS, eq. (|2Uj) . we get 

AH = -^H 3 . 

Therefore, from definition (|8d|) we obtain, 

SZ H = j^r e (3H 2 + 2Tt(H 2 )). 

With this value for SZh we can immediately find the one-loop (3 function for H, 

^=P = ~H + _L (3H 3 + 2Tt{H 2 )H) + 0(H 5 ). (30) 

/3, which does not depend on A at one loop, commutes with H. This need not be the case in more 
complicated theories where (3 can depend on other dimensionless matrices. 

With 8M and /3, from ()27|) and ()30|) . we obtain the evolution equation for M, 

H^r = -7» = {^HMH + \{M, H 2 }^j + 0(H*). (31) 

The matrix ~f m defined by this equation is unconventional in that it has mass dimension one. At 
one-loop level we can take M on the rhs to be diagonal (in fact, we can set M = -Mph) and set 
H to its tree-level value, neglecting its ~p dependence. The solution to (|3Tj) up to terms of 0(H 4 ) 
is then, 

M(/Z) = M(7Z ) - lm log ( P ) . (32) 
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We see that we can diagonalize the one- loop mass matrix at a given scale ~fi . If [M, H] = 0, then 
M also commutes with ■y m and we can diagonalize M at all scales. This is the case, of course, if 
M oc 1 or H oc 1. It is easy to check that the expression (|29[1 for M is a solution to (j31j) . 

We consider, finally, the anomalous dimensions of fields. From the expressions for SZ^ and SZ 
we find, 

Once again, the dependency of H on ~p can be neglected in the rhs of these equations. The evolution 
of if) with Ji is given by an hermitian matrix, which is consistent with U=l. 



3 Normal mixing matrix 

The treatment given in ^Hfor the case of a hermitian mixing matrix can be readily generalized to 
normal mixing matrices. Such a generalization is natural, since it comprises all mixing matrices 
that can be diagonalized by a unitary transformation of fields. Consider, for instance, a Yukawa 
interaction of the form, 

£ = -4{a + + ^ (i - M )V> + ? W ( A + ^ N\^A ~ V (</>, </> f ) , (34) 

where V contains cubic and quartic terms and [N, N'] = 0. This interaction induces 4>-(ft mixing, 
which must be taken into account when renormalizing the theory. Besides that extra complication, 
the renormalization can be carried out along the same lines as in where now, 

N = U^ N Z N NU N , (35) 

with U n unitary and Zn normal such that [Zn,N] = 0. Unlike the previous case, though, we 
cannot parametrize N in the form B^NB, since congruence transformations do not preserve 
normality. 

More generally, we must ask whether the normality of N is stable under renormalization. (The 
hermitianity of H in fj2]was obviously stable by unitarity.) That this is so can be seen by noticing 
that the Yukawa interaction in (|34|) is invariant under a (U(1)) N symmetry which, in a flavor 
basis in which N is diagonal, is given by ip a — > e~ taa ^ a , <fi — > 4>. This (U(1)) N symmetry will 
generally be broken by fermion mass terms. In MS the Yukawa coupling will remain normal 
after radiative corrections are taken into account, since the fermion wave-function renormalization 
matrix is diagonal in the interaction basis due to the (U(l)) symmetry (or, in this model, because 
it must be a polynomial in AT"). In OS finite asymmetric counterterms to the Yukawa coupling are 
needed, which render the renormalized Lagrangian asymmetric. An example of this well-known 
phenomenon (see, e.g., jlUl Ch.4] and refs. therein) is considered in the next section in connection 
with unitary mixing. 



4 Unitary mixing matrix 

The case of a unitary mixing matrix constitutes an important example of non-hermitian, normal 
mixing. Besides its obvious phenomenological relevance this case is interesting because 

mixing matrix unitarity imposes restrictive constraints on the form of counterterms, on top of those 
stemming from normality. 
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Furthermore, by considering unitary mixing we can extend our approach to interactions of the 
general form Vi^V^i with two different families of fermions. Independent unitary transformations 
of the fermion fields lead to a biunitary transformation of the mixing matrix V. In general normality 
is not preserved by biunitary transformations, so our previous results do not apply to this type of 
interactions unless V is unitary. 

In this case, the (U(l)) N flavor symmetry of the previous section is replaced by a larger SU(N) 
invariance, since we can choose the phases of fermion fields so that in the interaction basis the 
mixing matrix is just the identity. For concreteness we consider the simple Lagrangian, written in 
the interaction basis, 



£ = £ 7 + £ + £^ + C Y + C cm - V (36) 
\F, U F^ + \m 2 A,A» - 1 



£ 7 = — A F^F^ + -n^A^A* - - (d^f , £ = -<f>i(p + m 



2 



£em = QytftA^ + (-j; +j x + ft) A„, V=^ (fcp) 2 

3% = Q^ei (d^<j> - , # = Qje^^j . 

The model consists of two fermion families tpi 2 , each containing N flavors. These interact with a 
scalar field <p trough the Yukawa couplings in £y, which are arbitrarily chosen to be chiral, with 
P± = (1 ± 7s)/2. We choose g to be real and positive, since its phase can always be absorbed in 
(f>. All of these fields are charged, with Qs = Qi — Q2, and minimally coupled to a massive photon 
field A^. U(l) gauge invariance is broken explicitly by the photon mass and by the covariant 
gauge fixing terms. One- loop radiative corrections to the Yukawa interaction arise only from the 
couplings of if)j and 4> to the massive photon. 

Besides U(l) gauge invariance, £ has a global SU(N) flavor symmetry ipj — > e~ l ^ AaAXA iJ>j, 
broken by fermion mass differences, with associated current, 

2 2 

3 A = Y/^i'f^j, 9,Ja = i E^i (P+[Mj,X A ] + P-[M], \a}) fy. (37) 
3=1 3=1 

Here, (A = 1, . . . ,N 2 — 1) are the generators of SU(N) in the fundamental representation. £ 
has also broken global axial t/(l)s and SU(N)$ symmetries which will not be needed in the sequel. 

The renormalization of this model proceeds much in the same way as in ^2 We introduce 
renormalization constants Z^, SmZ for the <f) field and its mass. For the Yukawa coupling we 

write g = ^ e ^ 2 Z g g , and for the scalar self-coupling, A = //(A + 5X). By gauge invariance, only 
wave- function renormalization is needed for the gauge sector [SJ , 



A^ = Z l ' 2 A\ mL = ZX. & = £y£> e =fi t/2 Z- 1 / 2 



e. 



Henceforth, we work in a flavor basis such that at tree level mass matrices are diagonal. We 
can always find unitary matrices Vj, Wj, (j = 1,2), so that the mass matrices are written 
as Mj = VjM'jWj, with M'a diagonal, with positive entries ^2 §21.3]. The corresponding 



transformation of fermion fields is tpj = (WjP+ + VjP-)iJ)j. Writing the Lagrangian in terms of 
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primed quantities, and dropping the primes, at tree level we get, 
2 

= Y,4j(i9- M ) C y = O^iP+V^ + g^P-V^^, (38) 

3=1 

with V = V\W 2 the unitary mixing matrix. In this basis, the SU(N) global flavor symmetry is 
given by, 

%f> x — ► e ^ E ^' 4 ^Vi, ^2 — ► V ] e~^A a A\ A y^ 2 ( 39a ) 

f A = ^^Xa^i + ^V^X A V^ 2 , d„f A = i^jMi, AaJ^i + ty 2 [M 2 , V^X A V}^ 2 . (39b) 

Notice that we can parametrize this flavor symmetry, and its associated currents, in an infinite num- 
ber of different ways. We could have written, for example, ipj — ► (P + V 1 ^ 2 ^ exp(— i J2a 01 aXa)V 1 I 2 + 

P- V 1 ' 2 exp(-i Ea <*aXa) V 1/2t )Vi- 

Using the same argument as above, we can apply a finite unitary renormalization of fermion 
fields so that M\ and M 2 are hermitian in each order of perturbation theory. Assuming this has 
been done, we can use ()A.3j) to write, 

M j0 = ul jL (Mj + SM^U^r, (40) 

so that in 

(P+M j0 + P_M} ) = UmjiMj + 8M 3 )U mj , with (41) 
U mj = P + U mjR + P_ U mj L , U mj = 7°wLi7° , U rn3 u\ nj = 1 = u] nj U m3 . 

In OS these expressions are always valid since Mj are real diagonal to all orders, therefore hermi- 
tian. In MS we can use ()4()l) at one loop, because the tree-level Mj have been chosen hermitian. 
The resulting one-loop renormalized mass matrices will not be hermitian. As mentioned above, we 
can restore hermitianity at one loop, at a given renormalization scale, by means of a finite unitary 
renormalization of fermion fields. This is analogous to the choice of one- loop off-diagonal entries 
of the mass matrix in ^2.41 

The bare fermion fields are given in terms of renormalized ones by expressions similar to Q , 



^ j0 = UjZ)% (42) 

Z 1 / 2 = (P + Z 1 / 2 + P^Z 1 / 2 ), U J = (P + U, R + P-U, L ), UjU} = l. (43) 

The mixing matrix V introduced in (|38j) is renormalized according to (jA.lJ) . 

V = W^ZyVW, [Zy, V] = 0, W = e~ iSW . (44) 



Since Zy is unitary, we can write Zy = exp(-idZy), with 8Zy hermitian, [SZy, V] = 0. Other 
parametrizations for the mixing matrix counterterms are discussed in the next section. 

Substituting these expressions for bare quantities in C, we obtain its expression in terms of 
renormalized parameters and fields. The fermion and Yukawa Lagrangians, in particular, read 

+ C Y = Y^jZT* $ Z T^i ~ Y/^i^T^miM + SM j )U mj U j Z 1 / 2 ^ j 

j 3 ' (45) 

+ ^ 2 Z g zl /2 g^ 1 P + Z 1 1 /2 u\ L W^ZyVWU 2 RZl /2 tP 2 c) ) + H.c. 
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Just as in ^21 out of the unitary matrices 14 j, U m j, W that we have introduced, only the combina- 
tions U m jUj, WU\l and WU2R enter C. Clearly, we can always choose Uj = 1. In OS, however, 
it is convenient to set li m j = 1 instead, so that bare mass matrices are diagonal in the mass basis. 

We see from (j45f) that for the interaction term to retain its form, with a unitary mixing matrix, 
we need Z\l,Z2r oc 1. This is the case in MS , due to SU(N) flavor symmetry. In OS finite 
asymmetric counterterms to the Yukawa coupling are needed, so the form of the Lagrangian is not 
preserved. 

4.1 Other parametrizations for the mixing matrix 

The parametrization of counterterms given in (|44j) conforms to the general form for normal matrices 
given in (jA.lj) . Since V and V are both unitary, however, we can write the relation between them 
in other ways. Clearly, V = WV or V = VW are admissible since we can reach any unitary 
matrix in a neighborhood of V by varying W over a neighborhood of the identity in SU(N). 
Another usual way of writing the renormalization constants for V is V = W1VW2, with W\ t 2 
unitary 015- This parametrization is convenient from the calculational point of view. Here we 
point out that it is overspecified, it must satisfy constraint relations analogous to those considered 
in H2.ll as we show next. 

Given a map F : SU(N) -> SU(N) of the form F(V) = WiVW 2 , Wj(V) G SU(N), we can 
always find Wi 2 such that, 

F(V) = WtVW 2 , Wj = Wj{V) G SU(N), [W 2 Wi,V] = = [W{W 2 , F(V)}. (46) 

To see this we notice that, given Wip, we can use Lemma U of appendix E] to write F(V) as in 
(jSJ) (or (|A.ljl ). With the same notation as in eq. (JMJ), setting W\ = W^Z V and W 2 = W we 
get (gHJ). 

4.2 On-shell scheme 

We consider only the case of regular mass matrices Mj. The extension to degenerate Mj can 
be carried out as in Renormalization conditions in this scheme break flavor symmetry, so 

counterterms are not symmetric, their finite parts being tuned so the field basis is such that Mj 
are diagonal to all orders. Setting U m j = 1 in (|4*T1) we have Mj = Mj + 5Mj, with SMj diagonal. 
To one-loop order we then have, 

j j 

- Y$3 M rfj - Y^j (P+ AM J + >' AA/ ,) *j 

j j 

C Y = ^ l2 g$ x P+ (V + AV + VAg + AT) V> 2 <£ + H.c. ^ 
C em = Ql^e 2 (1 + Re(5Z^)) <fi <f>ApA» - // /2 j£ (1 + Re(5Z^)) A„ 

+ Yl M e/2 i; ^ + E Qj^ 2 ^^ (P-Sz.l + p+sZjr) ^ 3 a^ 
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where the counterterms are defined as, 

AMj = SMj + idU jL Mj - iMjdUjR + ^SZ jL Mj + ^MjSZ jR 

AV = i[SW, V] - iSZ v V + iSU 1L V - iV6U 2 R (48) 

Ag = 5Z g + ±5Z (j> ; AT = Uz lL V + ^VSZ 2R . 

The flavor-asymmetric counterterms must be finite, divergent terms being flavor symmetric. In 
particular, the divergent parts of wave-function renormalization (SZj^^)^ oc 1. The phase of 
8Z<p appears only in Ag, so it can be adjusted to keep Ag real. We defined AV so that it is 
perturbatively unitary. In a flavor-symmetric, mass-independent scheme such as MS , fermion 
wave-function renormalization constants are flavor scalars that can be absorbed in Ag, no other 
counterterms to the Yukawa coupling being needed. In particular, there are no flavor-breaking 
dimension 4 operators in C In OS, however, we also need AT, whose finite part can be viewed as 
either breaking the unitarity of V or the scalar nature of g. 

For the computation of Feynman graphs involving fermion loops with 75 vertices we use 't Hooft 
and Veltman's prescription, with 75 anticommuting with 7^ for fi = 0, . . . , 3 and commuting other- 
wise. We henceforth set = for the sake of simplicity. For the scalar field self-energy, with OS 
renormalization conditions, we find, 

IF> 2 ) = ^(p 2 ) - fy(m$) - (p 2 - m 2 )n'4m 2 ) (49) 
g 2 N (2 



e 




5m l = -JKZz ( 7 ( J2 m % ~ ~k m \ I + -T m l ~ ^Zrfat 1 + a o( m ja)) I + 



tt<f>(p 2 ) = -^2^Z V ^V\ a {p 2 - m ia -m 2 2b )b (p 2 ,m 2 la ,ml b ). 

a,b 

The fermion two-point function for each family can be expressed in terms of form factors as 
Tjip 2 ) =p! -Mj-TL^p 2 ), 

Ujip) = (j/£jV) + AAft)P_ + (j*Ef (p 2 ) + AMj)P + + Sf (p 2 ). 
The form factors for the first family are given by, 



^labiP 2 ) = - Sz ub ~ 77TT 6 ab ~ Qla 2 ^^5 ab + ti( ab (p 2 ) + Q\^ m (p 2 ,m; 2 ,m\ a )5 ab 
^ ab (p 2 ) = -bZ* ab - Qle 2l ^ ab + Qi^ m (p 2 , m 2 ,m 2 a )5 ab (51) 



^fab(P 2 ) = ™iaQie 2 ^2~<U + m la Qln^ m (p 2 , m 2 , m 2 la )5 ab . 



It is not difficult to check that the dependence on £ in the scalar and fermion propagators satisfies 
U(l) e . m . Ward identities to one-loop level |11[ §18.7]. Notice the flavor dependence of e.m. correc- 
tions in (|51|) . This flavor structure is also apparent in the counterterms to the e.m. current in (|47|) . 
In equation (|5Tj) we have introduced the quantities, 

VU(P 2 ) = ^Yv ac vib.(p 2 ,m 2 ,m 2 c ) (52a) 

c 
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&^(j> 2 ,m*ml) 



16vr 2 



l + 26_(p 2 ,m 2 ml) -b (p 2 ,m 2 m \) +£b (p 2 ,£ml, w 2 ) (52b) 



p — m c 



^em(P 2 > m 7' m ? 



16tH 



(2 + 35 (p 2 , m 2 , m 2 c ) + £,b (p 2 ,£m 2 , m 2 c )) 



(52c) 
(52d) 



The corresponding expressions for the second family are obtained by changing the family index 
1— >2 and L<->R, with, 



16vr 2 



The renormalization conditions are 



mA a { P 2 ) + m ja Xf aa (p 2 ) + 2Y, s jaa {p 2 



(52e) 



0, (53a) 



for diagonal functions and 

r-R.L 



m^f (m l6 ) + AAf*f + Sf a6 (mf b ) = 0, m ia E^f (m? a ) + AMjf + Sf a6 (mf a ) = 0, (53b) 

for both diagonal and off-diagonal ones [1]. Here we have used the notation AM R ab = AM jab and 
AAf^j = AMt 6 for brevity. The renormalization constants are then, in the flavor-diagonal case, 



5Z R 



laa 



szL 



AM 



laa 



Ql e £ . n 2 n V ^ 2 2 2 \ , 2 o^" ^ 2 \ 

— + Qi^em {m la , m m la ) + m la il laa (m la ) 

+ 2 rnl a Qln^(ml a ,m 2 ,ml a ) + 2m 2 la Q 2 1 VL s e ' m {m 2 la , m 2 , m 2 la ) 

rtfriS I 2 2 2 \ r&c^V I 2 2 2 \ , jfzJ? 



laa • 



In the off-diagonal case e.m. contributions vanish and we obtain the simpler expressions, 



~^la m lb 



m la ~ m 
,2 nY /_2 



AMi a6 



™l<Xa&(™la) " m lb^iab ( m lb) 



m 



la 



m 



lb 



m la^Z? h . 



(54a) 

(54b) 
(54c) 

(54d) 

(54e) 
(54f) 



We also obtain AM^ ab = m\ bZ R ah , which is consistent with the above definitions. Notice that 
off-diagonal counterterms are finite, as required by flavor symmetry. Substituting these results 
back into (|51[) we obtain the renormalized fermion self-energy for the first family. 

In order to compute SU^' R we proceed as in It is convenient to split AM into its hermi- 
tian and anti-hermitian parts and, using (|48[). write an equation analogous to (|16|) for each part. 
Proceeding in that way we get, 



SM lab = 6 ab ^ (5Z* a - bZ[ aa ) 



(55a) 
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2vfi\ a -vfi\ h rn\ a -m\ b 



su l _ ■ m la m lb R i mj a + m 2 lb L 
m ia ~ m( b 2 mf a - mf b 



where in the last two lines a ^ b. Diagonal elements of SU 1 ' are not determined by renormaliza- 
tion conditions, and are set to vanish. Renormalization constants for the second family are given 
by similar expressions, after changing the family index and exchanging L and R labels. 

We denote the three-point 1-PI function corresponding to the Yukawa vertex as T^{p\,p 2 ) 

l 2b' 

3/*7 /4 + 1 >!/< ^ ^ J - :{,<;-' 

Only Tt receives divergent contributions at one loop, 



External momenta are assumed to be on their mass shell, p\ = m\ a , p\ = m 2 b , and = (pi+p 2 ) 2 — 
ml. T 3 can be decomposed in form factors as T 3 = T+P+ + T~P_ + tY ^ + Tg 4 7^75 + Y\ . 



T t,ab = 9(V + VAg + AV + AT) ab + ge 2 Q 2 V ab l(p u p 2 ; m la ,m 2b ), (56) 

where the loop integral I is defined in appendix [D] We separate in Ar the contributions coming 
from flavor-diagonal fermion wave-function renormalization constants, which are divergent and 
gauge-dependent, from the finite, ^-independent off-diagonal ones. Thus, 

rj ab = g(V + VAg + AV) ab + g(AT) + gV ab Q ab (57a) 

V / ab 

( Ar ) ab = \ E 6Z lacVcb + \Y1 V ^ 6Z 2cb (57b) 
c^a cj^b 

Gab = \^>Z\ aa + -8Z2 bb + e 2 Q 2 2(pi,p 2 ;m la ,m 2b ). (57c) 
An explicit expression for Q ab is as follows, 

n 9 2 . 3Q 2 e 2 1 / Y / 2 \ 1 2 r>y t 2 \ 1 / 2 2 \ 

= ~\^{T e + "8^7 2 V laa ^ + m la fi laa( m la) + i m la ~» ™2b) 

+ ^2 (-^ A ( m ia>^)-^A(ml 6 ,m 2 ) + A e (m? a ,m 2 ) + A e (m| b ,m 2 ^ ! of>) 



where, 



+ 3Q e p 1 -p 2 C {p (j) ,p 1 ;m 1 ,m la ,m 2b ) + Q e p x ■ p 2 £Co{p (j) ,p 1 ]£ > m 1 ,m la ,m 2b ), 

A(m 2 a , w 2 ) = 1 + fri(m 2 a , m 2 , m 2 a ) + 2b (m 2 la ,m 2 , m 2 la ) + 4m 2 la b' (m 2 la , m 2 , m\ a ) 
+ Am 2 la b' 1 {m 2 la ,m 2 ,m 2 la ), 
?n 2 

A^(m\ a ,m 2 ) = — ^ (6i(m 2 a ,m 2 ,m 2 a ) - 6i(m 2 a , £m 2 , m?J) , 



(59) 



and p\ ■ p 2 = 1/2(771? — m 2 a — rn^ b ). Some comments regarding Q ab are in order. The ultraviolet 
divergent terms in (|58|) are ^-independent, as expected since they must be cancelled by Ag. All 
the remaining gauge dependence in Q ab is contained in A^ and C$. As is well-known, and is clear 
from its definition in appendix El the triangle integral Co with on-shell external momenta diverges 
logarithmically as m 2 — ► 0. So does A^ as well, since we can write, 

A ? (m 2 , m 2 ) = i (b (m 2 ,m 2 , m 2 ) - a (m 2 )) - - (bo(m 2 , £m 2 , m 2 ) - a (£m 2 )) , 
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with ^ 

bo(m 2 ,m 2 m 2 ) — ao(m 2 ) = / (fx log ( 1 — x + ^ttX 2 ) + 1. 
7 7 Jo V «7 / 

These infrared divergences must be cancelled in the computation of transition rates by the con- 
tribution from real-photon emission diagrams. We conclude that the infrared divergent terms 
must appear explicitly in the amplitude, and should not be absorbed in the finite part of coun- 
terterms. In fact, the requirements of infrared and ultraviolet finiteness, that g and g should be 
flavor-independent and V and Vo unitary, and of gauge invariance, lead us to set, 

a 2 30 2 e 2 

A 5 = -(<?afc)div = y^T € ~ ^27> and AF = 0. (60) 

These equations determine the renormalization constants 5Z g , SW and SZy through (|48j) and 
the above results for wave-function renormalization constants. Our choice is not unique, though, 
since the above requirements still allow the finite part of 5Z g to be redefined by adding gauge- and 
flavor-independent arbitrary constants. We could, for instance, define g through the total <j) decay 
width, or some other inclusive process. We will not pursue such phenomenological analysis further 
here. 



4.3 MS scheme 

Renormalization in MS scheme turns out, as expected, to be much simpler than in OS. The 
renormalized Lagrangian at one-loop level takes essentially the same form as in (|47|). (|48[). The 
main differences being that now SU^ ,R = 0, SZj ,R = 5Z^ ,R 1, and SU^f 7^ 0- The mass matrices 
Mj are diagonal at tree level but not necessarily diagonal, or even hermitian, at one loop. We are 
then led to separate diagonal elements (which are real, given our choice of tree-level flavor basis) 
from off-diagonal ones (which are 0(g 2 ) or 0(e 2 )), as discussed in ^2.41 

We define Ag = 5Z g + 1/2SZ^ + 1/2SZ[ + 1/25Z?, and set AT = 0, since obviously no 
asymmetric counterterms are needed in this scheme. Furthermore, SW = for the same reasons 
as in H2.41 Renormalization constants can then be read off the corresponding expressions in H4.21 
In particular, 

h7 L Q 2 e 2 C 9 2 , 7 r R? 9 2 (N + 1) 3Q 2 e 2 

Divergent contributions from the triangle diagram correcting the Yukawa vertex, which are flavor 
symmetric, have been absorbed in SZ g . There are no further infinities to cancel, so that SZy = 0, 
i.e., V is not renormalized. 

We can easily derive renormalization group equations for renormalized parameters and Green's 
functions. However, V obviously does not run at one loop in this model. 



5 Concluding Remarks 

Unitary transformations play a distinguished role in quantum theories, as is well known. In QFT, 
unitary transformations in some internal "flavor" space preserve the normalization of kinetic terms 
and currents, in particular electromagnetic ones. It is then natural to consider the class of mixing 
matrices that can be diagonalized by a unitary transformation, namely, normal mixing matrices, 
and their renormalization properties. This is the subject of the foregoing. 
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We point out in £J21 and |JH that the usual parametrizations for mixing matrix counterterms 
HI El E] are over specified. This fact, of course, does not make them less useful in any way, but 
had not been pointed out in the previous literature. We explicitly exhibit the constraints those 
parametrizations must satisfy, and their solution. Furthermore, the minimal parametrization given 
in appendix E] can be applied also to the case of normal mixing matrices, which we do in 

In principle, the general case of non-singular, not necessarily normal mixing matrices, although 
not discussed in this paper, can also be studied in the framework proposed here for normal matrices. 
This is so because any complex matrix A can be decomposed as A = N + T, with N normal 
and T nilpotent such that they can be simultaneously brought to diagonal and strictly triangular 
form, respectively, by a unitary transformation. 

Writing field-strength renormalization matrices in polar components is useful in OS scheme, as 
discussed in §2 and §4. At one-loop level this is the same as decomposing those matrices into their 
hermitian and anti-hermitian parts, as done in pQ. As remarked in §2, the unitary and hermitian 
components play different roles in the theory. Furthermore, the contribution to higher orders in 
perturbation theory from powers of the one-loop SU is apparent in this way due to the exponential 
form of the unitary matrix U = exp(i8U). 
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A Mappings of normal matrices 

A matrix N £ C NxN is called normal if [JV.jyt] = [T3]. It can be shown ^3] that a complex 
matrix can be diagonalized by a unitary transformation if and only if it is normal. We denote by 
n(N) the set of N x N normal matrices. The Lie algebra u(N) of hermitian matrices, the group 
U(N) of unitary matrices, and the Abelian algebra of diagonal matrices are all contained in n(N). 

Lemma 1. Let F : n(N) — > n(N) be a mapping of normal matrices such that rank(F(X)) < 
rank(X) for every X in n(N). Then, we can write, 

F(X) = U ] ZXU, with [Z, X}=0, Z = Z(X) G n(N), U = U(X) G U(N). (A.l) 

Proof. We first consider the case of X non-singular. In fact, in this case there is an open neighbor- 
hood of X in n(N) where all matrices are non-singular. We can apply the parametrization (|A.1|) to 
every matrix in that neighborhood. Using the normality of F(X) and X , we can write (dropping 
the argument X for brevity) F = V^F'V and X = W^X'W, with V, W unitary and F' , X' 
diagonal. Define Z' diagonal by F' = Z'X'. Then U = W ] V G U(N) and Z = W ] Z'W G n(N) 
satisfy (|A.1|) . 

Assume now X has a null eigenvalue with multiplicity r. We can choose W so that X' = 
diag(0, . . . , 0, x r +i, . . . , x n ) with Xj ^ 0. Similarly, we can choose V so that F' = diag(0, . . . , 0, f q +i, 
. . . , f n ), with fj 7^ and q > r. We then define Z' = diag(l, . . . , 1, f r +i/x r +i, ■ ■ ■ , fn/x n ) and pro- 
ceed as above. 

A transformation F satisfying the rank hypothesis of the Lemma may be called "multiplicative," 
in the sense that F(0) = 0. The same argument as in the multiplicative case, with obvious changes, 
proves the following. 

Lemma 2. Let F : n(N) — > n(N) be a mapping of normal matrices. Then, for any X G n(N) we 
can write, 

F(X) = U\X + N)U, with [A/", X] = 0, AT" = N(X) G n(N), U = U(X) G U(N). (A.2) 

We remark that both (jA.l|) and (jA.2|) hold for mappings of hermitian matrices, F : u(N) — > 
u(N), with Z and N hermitian, and ()A.1|) also holds for mappings of unitary matrices, F : 
U(N) -> U(N), with Z unitary. 

Consider, finally, a general matrix function of a hermitian matrix F : u(N) — ► C NxN . Using 
the polar decomposition F(X) = R(X)U f(X) with U f unitary and R hermitian, and applying 
eqs. HA.1|) and QA.2J) to R, we obtain the parametrizations, 

F(X) = UZXV, U = U(X) and V = V(X) G U{N), Z = Z(X) G u(N), [Z, X] = 0. 

(A.3a) 

F(X) = U(X + SX)V, SX = SX(X) G u(N), [SX,X] = 0. (A.3b) 
Substituting (|A.3|) into the hermitian matrices FF^ and F' F, we recover ()A.1|) and (|A.2|) . 
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B Perturbative factorization of SU(N) matrices 



In ^21 we make use of the fact that a unitary matrix U which is close to the identity can always 
be factorized as U = exp(iSU') exp(i5U), where SU' is hermitian and diagonal and SU is her- 
mitian and has zeros on the diagonal. A sketch of the proof follows (see |14j for a globally valid 
factorization) . 

We decompose the Lie algebra of hermitian matrices as u(N) = a © b, where a is the Cartan 
subalgebra of diagonal matrices of u(N) and b its complementary subspace. Let U = exp(ieif), 
with e small and H G u(N). We want to show that we can always find A = A(e) G a and 
B = B(e) G b such that 

N (■ \n 

^eA e ieB = ^2^1L H n + 0(e N+1 ). (B.l) 
n=0 n ' 

To this end, we write A = z~2n=o nl^-n with A n G a, and analogously for B, and consider the 
equation exp(ieH) = exp(ieA) exp(ieS) which, using the Campbell-Baker-Hausdorff formula (see 
[HI E] an d references therein) can be written as, 

H = (A + B) + B) + ([[A, B],B] + [[B, A], A]) H (B.2) 

Expanding the rhs in powers of e, we are led to a set of recursive equations, 
A + B = H 

A 1 + B 1 = -[B ,A ] (B.3) 

A 2 + B 2 = [BuAq] + [B , A!] - h[Ao,B ],B ] - h[B ,A ],A Q ], etc. 

fa fa 

which can be solved iteratively up to the desired order. (Notice that in ()B.3|) we have used 
the Abelianity of a to eliminate some commutators.) Thus, Aq = diag(f/n, . . . , H^n) and Bq = 
H — Ao, A\ = and B\ = 1/2[Bq, Aq], and so on. More generally, the equation for the coefficients 
of order iV + 1 obtained from (jij.2|) is of the form, 

N 



... (A N+1 + B N+1 ) + - -^- } [A m ,B n ] 
+ 1)! 2 ^-^ „ mini 



(N + 



m,n=0 
m+n=N 



1 iV_1 1 

+ 7^ E 1 - r - n ([[A k ,B m },B n }-[[A k ,B m ],A n ]) + --- 
12 , „ k\m\n\ 

k,m,n=Q 
k+m+n=N-l 

where the ellipsis refers to higher-order commutators. Terms with if -order commutators involve 
Aj, Bj with j = 0, . . . , N+ 1 — K, which are already known. Thus, by projecting (A^v+i + Bjv+i) 
as given by this equation over a and b we find A^v+i and Sj\r+l- 

To summarize, let D be a neighborhood of the identity in U(N), determining a corresponding 
neighborhood d of in u(N), where the CBH formula holds. Then, the elements of the form e A e B 
with A £ a and BEd are dense in D. 
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C Kernel and image of the adjoint map of a normal matrix 

Consider a fixed N x N normal matrix N. Then, any matrix A 6 C JVxJV can be written as 

A = B + [N,D] with [N,B]=0. (C.l) 

Equations Ijl6|) and (|21|) are of this type. Given A and AT, (jC.lj) can always be solved for £? and 
Z), as we now show (see also [HI §3.3]). 

Since N is normal, it can be diagonalized by a unitary matrix. Let Ai,...,A r (r < N) be 
its eigenvalues and dj(j = 1, . . . , r) their multiplicities. The associated orthonormal basis of 



of eigenvectors of AT" is denoted by {|Aj, ay)} (j = 1, 
Eto AilAi.aXAi.al, and for any A G C NxN , 



, r, Qj = 1, • • • , dj). Then we have N 



A= ^\\ h a){\ h l3\. 



(C.2) 



We define the adjoint map associated to N as adjy : C NxN -> C 7Vx7V , adjv(A) = [N, A]. Then 
the set of orthogonal projectors {|Aj, ati)(\j, f3j\} is a basis of C NxN of eigenvectors of ad^, since, 
[N, \Xi, cti)(\j, = (A, — Aj)|Aj, ai){Xj, /3j\. Therefore, we can rewrite (|C.2|) as 



i,a,/3 



A™\\ i ,a)(\ i ,P\ + 



1 



"'A, A 



|Ai,a)(Aj,/?| 



(C.3) 



where the first term on the rhs obviously commutes with N. This is the decomposition (jC.lj) . 

Whereas (jC.3|l provides an explicit solution to (|C.1|> . a slightly broader point of view can 
be adopted. Consider the inner product (A, B) = Tr(A^B) in C NxN . By the cyclic property 
of the trace, (A,&d N (B)) = (A,[N,B]) = Tr (A* [N,B]) = Tr{[N*, A]*B) = ([N*,A],B) = 
(ad^t (A) , B) , so (adjv)^ = adjyf With this, using the Jacobi identity and the normality of AT, it 
is immediate that [(adjy) , ad^v] = and thus ad;v is a normal transformation of C NxN . Therefore 
the spectral theorem ^3] holds for adAr. In particular, C NxN = Ker(adAr) © Im(adjv)- This 
orthogonal decomposition shows that solutions to (|C.1|) exist and are unique up to addition to D 
of a matrix commuting with N . 



D Loop integrals 

In this appendix we give a list of loop integrals used in the foregoing. More complete calcula- 
tions can be found, e.g., in [1J El UB El- Divergent integrals are separated in a dimensional- 
regularization pole term and a finite remainder. ~p = fiV 47re~ 7B . 

A (m 2 ) 
a (m 2 ) 
A^(m 2 ) 



ifj, 

(2tt) £ 
/ m 2 

fog =J 



Ifl 

(2vr) c 



1 

m 2 + ie 



p. _ m 2 _|_ £ £ 



m 

8^ 



+ 



m 
16tv 



;a (m 2 ) 
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a 2 (m 2 ) = a (m 2 ) - - 
B (p^mlml) = ld% ' 



(2ir) d 7 (£ 2 - m\ + ie) ((£ + p) 2 - m 2 , + ie) 
' +^^o(/,m 2 ,m 2 ) 



8tT 2 € 167T 2 



&o(p 2 ,m 2 ,m 2 ) 



f~^~ ( 777? 777? Jp 1 \ 

/ cfa log I (1 — ac) + — x(l — x) ^ — ie I 
V ^ M M / 



(2vr) d J (£2 - m \ + ie) ((£ + p) 2 - m\ + ie) 



P M . / 2 2 2^ 

16^-16^ 6l(P ' mi ' m2) 



6i(p 2 ,m 2 ,m 2 ) 



which can also be written, 



/ dxx log I (1 — x) + — x{\ — x)=2 — ie J 



2 2\ 

,m5) 



B^(p^,ml,ml) =p^Bi(p 2 ,m\ 
p 2 B 1 (p 2 ,m\,ml) = - (A {ml) - A {ml) - (p 2 + m\ - ml)B (p 2 , m\, m\)) 

We also use the combination b-(p 2 , m\ , to 2 ,) = bo(p 2 , m 2 , m|) — bi(p 2 , m\ , m 2 ). 

C (pi,P2,ml,ml,ml) = -^—r [ d d £ 



(27r) d 7 (£ 2 -m\ + ie) ((£ + p 2 ) 2 -m 2 + ie) ((£ - pi + p 2 ) 2 - mf + ie) 
Co is ultraviolet finite. In 82,21 we use the following triangle integrals, 

1 i / +p , 2 + m 2 £f- $1+^2 +m 3 



Hi(pi,p2]m 1 ,m 2 ,'m 3 ) = , n / d £—^ 5 — — — r~ 5 - — ^ * — — 

v 1 2 d/ (2vr) a! J £ 2 -mj + ie(£ + p 2 ) 2 -ml + ie(£-p 1 +p 2 ) 2 -ml + ze 

Tr(Hi(p l ,p2;mj,m 2 ] ,ml)) = Tr(l) ^-^^ + h 1 (p 1 ,p 2 ;ml,m 2 ] ,ml) 

h^p^mlrnlm 2 ) = \b Q (( Pl - p 2 ) 2 ; m 2 , m 2 ) + \b Q (p 2 ;m 2 ,m 2 ) 

+ \ (K + m 3 ) 2 -p 2 ) Co (pi,ja;m?,m2,m§) 
h 2 (pi,P2;m 2 ,ml) = Tr(l)m 2 C (pi,p2;m 2 ,,ml,m 2 ) 

with Tr(l) =4. Finally, in 84.21 we define, 

Tfpi ^m 2 m 2 m 2 )- 1 ^- I d d £A (£) T " i^l± tl + ^) P ^~ V 2 + ^ b )Y) 
l( Pl ,p 2 , m la ,m 2b , m 7 ) - 2 ^ y d £A^t) {{£ + ^ + ^ ((/ _ ^ _ ^ + ^ 
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= m ~ 8^ ~ 3^2 ( 36 0(Pl> m ?> m lo) + 3& o(P2> m r m 26) +C&o(P?,^7^1a) 

+£&o(P2>£ m 7> m l>)) +Pi -P2 (3C o (p0,pi;m^,mf o ,m2 6 ) + £Co(?ty,Pi; £™ 7 , ™|>)) 
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Erratum 



Mixing-matrix renormalization revisited 
Eur. Phys. J. C 20, 239-252 (2001) 
q ■ Antonio O. Bouzas 

o 

(N 

In a recent article [1] we discussed the renormalization of normal coupling matrices. There are 
two erroneous remarks in [1] that we would like to correct in this note. 

In section 4 of [1] we consider the case of a unitary mixing matrix in the context of a model 
involving two families of N Dirac fermion fields each, with N arbitrary, coupled to a scalar and a 
■ massive vector field. If the renormalized mass matrices for fermions Mj, j = 1,2, are chosen in 

such a way that they are hermitian in each order of perturbation theory, we can write the relation 
among bare and renormalized mass matrices in the form (see (40) in [1]), 

M j0 = ul jL (Mj + SMj)U mjR , (1) 

with 

Ph! [M j ,SM j ]=0, (2) 

as shown in appendix A of [1]. This is always the case in OS scheme, in which Mi 2 are required 
to be real diagonal at any perturbative order. 

It is asserted in [1] that (1) and (2) hold also in MS and related schemes at one loop if the 
tree-level flavor bases have been chosen so that M\ t 2 are hermitian at tree level. That statement 
is incorrect. It is not difficult to show that, given the polar decomposition Mj = RjVj with Rj 
hermitian and positive and Vj unitary, the relation (1) holds with SMj = SRjVj where SRj is 
hermitian and [Rj,SRj] = 0, that is, 

M jM j.SMjSm) = = [MtMj^MpiMj], j = 1, 2, (3) 

instead of (2). The substitution of (2) by (3) in the analysis of the model in MS scheme does not 
affect the other results given in [1] in any way. 

Also in section 4 of [1] it is stated that due to gauge invariance only a wave-function renormal- 
ization counterterm is needed in the vector field sector of the model. This is true in MS scheme but, 
of course, other, finite counterterms are needed in OS scheme. The vector boson mass and gauge 
parameter renormalization, however, are not needed in the analysis presented in that section. 
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